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Abstract —We propose and investigate a compressive architec¬ 
ture for estimation and tracking of sparse spatial channels in 
millimeter (mm) wave picocellular networks. The base stations 
are equipped with antenna arrays with a large number of 
elements (which can lit within compact form factors becanse 
of the small carrier wavelength) and employ radio frequency 
(RF) beamforming, so that standard least squares adaptation 
techniques (which require access to individual antenna elements) 
are not applicable. We focus on the downlink, and show 
that “compressive beacons,” transmitted using pseudorandom 
phase settings at the base station array, and compressively 
processed using pseudorandom phase settings at the mobile 
array, provide information sufficient for accurate estimation of 
the two-dimensional (2D) spatial frequencies associated with the 
directions of departnre of the dominant rays from the base 
station, and the associated complex gains. This compressive 
approach is compatible with coarse phase-only control, and is 
based on a near-optimal sequential algorithm for frequency 
estimation which can exploit the geometric continuity of the 
channel across snccesslve beaconing Intervals to reduce the 
overhead to less than 1% even for very large (32 x 32) arrays. 
Compressive beaconing is essentially omnidirectional, and hence 
does not enjoy the SNR and spatial rense benefits of beamformlng 
obtained during data transmission. We therefore discuss system 
level design considerations for ensuring that the beacon SNR 
is sufficient for accurate channel estimation, and that inter-cell 
beacon Interference is controlled by an appropriate reuse scheme. 

Index Terms —Compressive, mm wave, 60 GHz, picocells, RF 
beamforming 

I. Introduction 

The explosive growth in demand for wireless mobile data, 
driven by the proliferation of ever more sophisticated hand¬ 
helds creating and consuming rich multimedia, calls for orders 
of magnitude increases in the capacity of cellular data net¬ 
works il]. Millimeter wave communication from picocellular 
base stations to mobile devices is a particularly promising 
approach for meeting this challenge because of two reasons. 
First, there are huge amounts of available spectrum, enabling 
channel bandwidths of the order of GHz, 1-2 orders of mag¬ 
nitude higher than those in existing systems at lower carrier 
frequencies. Indeed, channel bandwidths could potentially 
increase even further with advances in transceiver technology 
such as bandwidth/power consumption/linearity tradeoffs for 
ultra high-speed analog electronics, and speed/precision/power 
consumption tradeoffs for analog-to-digital converters. Sec¬ 
ond, the small carrier wavelength enables the realization of 
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highly directive steerable arrays, with a large number of 
antenna elements, in compact form factors, thus significantly 
enhancing spatial reuse. In this paper, we address fundamental 
signal processing challenges associated with channel estima¬ 
tion and tracking for such large arrays, placed within the 
context of system design for a mm wave picocellular network. 

While the signal processing and system design concepts 
presented here are of rather general applicability, our numer¬ 
ical results are for a particular setting that we feel has great 
promise, as also discussed in some of our prior publications 
El, El- We propose to employ the 60 GHz unlicensed band 
for base station to mobile communication in outdoor picocells. 
The base stations can be deployed on lampposts, rooftops or 
ledges, and have multiple “faces,” with each face containing 
one or more antenna arrays. An example deployment on 
lampposts in a zig-zag configuration (successive base stations 
on opposite sides of the street) along an urban canyon is 
shown in Figure [T] At the base station, we consider very 
large 32 x 32 arrays (such 1000-element arrays are still only 
palm-sized at a carrier wavelength of 5 mm) targeting the long 
term, as well as “moderately sized” 8x8 arrays (which can fit 
within an area of about half a square inch) which are currently 
realizable. Note that 16-element arrays were reported several 
years ago a, and are already deployed in existing 60 GHz 
products, while 32-element arrays have been prototyped 0. 
We assume that mobile devices are equipped with smaller 4x4 
antenna arrays. We focus on 60 GHz in order to leverage the 
significant advances that have occurred over the past few years 
targeting indoor wireless networks: once 60 GHz transceivers 
are embedded in mobile devices, one could use them to extend 
coverage to outdoor picocells, albeit with different approaches 
to medium access than in standard indoor wireless networking 
protocols. We focus on downlink 60 GHz communication, with 
the goal of enabling base station arrays to perform transmit 
beamforming towards mobile devices, despite the challenges 
posed by mobility and blockage (which occurs more easily 
at smaller wavelengths). We do not count on reciprocity. The 
uplink could be a standard LTE or WiFi link at lower carrier 
frequencies, used both for uplink data (not modeled here) 
and feedback for enabling spatial channel estimation at the 
transmitter. 

Multiple antenna systems at lower carrier frequencies have 
a relatively small number of elements, each with its own 
radio frequency (RF) chain. This provides control of the 
individual baseband signals associated with each element, en¬ 
abling sophisticated adaptation, including frequency-selective 
spatiotemporal processing (e.g., per subcarrier beamforming 
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Fig. 1; Picocellular network deployed along an urban canyon 
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Fig. 2; Channel sounding scheme: The transmitter repeats the 
M transmit beacons L times so that receivers can measure the 
channel Vij between every virtual transmit-receive pair 


in OFDM systems). This approach does not scale when we 
have a large number of antenna elements packed into a tiny 
form factor. Instead, we consider RF beamforming, in which 
a common baseband signal is routed to/from the antenna 
elements, and we can only control the amplitude and phase for 
each element. Indeed, we go even further, assuming that the 
amplitude for each element is fixed, and that we can only apply 
coarse four-phase control for each element. Standard least 
squares array adaptation and channel estimation techniques, 
which require access to the baseband signals associated with 
each element, do not apply in this setting. Instead, we consider 
here a compressive approach which exploits the sparsity of the 
mm wave channel, so that relatively few measurements are 
required for channel estimation despite the large number of 
array elements. 

Contributions: Our contributions are summarized as follows: 
Architecture: We propose a novel architecture in which base 
stations send out compressive beacons, with a different set 
of pseudorandom phases used to transmit each beacon. Each 
mobile measures the complex gains associated with each 
beacon compressively, using pseudorandom control of the 
phases of its receive array. The scheme, described in more 
detail later, is depicted in Figure |2l the base station sends M 
beacons, repeated L times, which permits the mobile to use 
L different settings of its own array to measure the associated 
complex gains. Each mobile feeds back appropriately chosen 


functions of its measurements to the base stations on the uplink 
(which can be a standard LTE link). Each base station use this 
information to estimate and track the dominant paths to each 
mobile that it receives feedback from. 

Algorithms: Eor the regular two-dimensional (2D) arrays con¬ 
sidered here, directions of arrivaFdeparture map to 2D spatial 
frequencies. The base station estimates the spatial frequencies 
to each mobile using a simple sequential algorithm, shown 
to be near-optimal (in terms of approaching the Cramer-Rao 
Bound) in related publications. The algorithm exploits the 
geometric continuity of the channel across successive beacon¬ 
ing intervals to reduce the required number of compressive 
measurements. 

System Design: While we do not provide a complete system 
design centered around our compressive strategy, we do pro¬ 
vide preliminary results addressing some of the most important 
issues. We show that the overhead incurred by our beaconing 
scheme is very small (less than 1%). Eurthermore, while com¬ 
pressive beacons are essentially omnidirectional (in contrast 
to the highly directive beams employed for communication), 
the link budget suffices for accurate channel estimation, and a 
simple beacon reuse strategy suffices to control inter-beacon 
interference across picocells. 

II. Background and Related Work 

There is a significant body of recent literature on the poten¬ 
tial for mm wave communication for next generation mobile 
cellular networks 0, Q, HI, a. Go), a. Eor these outdoor 
settings, most of the attention has focused on bands other than 
60 GHz; for example, Go) studies blockage probabilities and 
achievable throughput based on measurements in the 28, 38, 
71-76 and 81-86 GHz bands. However, as pointed out in our 
prior work a, 60 GHz has immense potential at the relatively 
short ranges of interest in urban picocells, with the oxygen 
absorption characteristic of this band only having a marginal 
impact on the link budget (e.g., only 1.6 dB at 100m). While 
studies such as a, a, a, m focus on coarse-grained 
channel statistics (and their implications on system design and 
performance), our focus here is on signal processing for fine¬ 
grained channel estimation. 

We assume that the channel is well described by a relatively 
small number of discrete rays with delays and angles of 
arrival/departure taking values in a continuum. The key contri¬ 
bution of this paper is to provide a super-resolution framework 
for estimating and tracking these rays, with model-based 
estimation allowing us to go beyond (spatial) bandwidth based 
resolution limits. To the best of our knowledge, the existing lit¬ 
erature on mm wave measurements does not attempt to super¬ 
resolve channels in this fashion, hence we do not know, for 
example, whether the continuous power-delay profiles reported 
in m are consistent with a parsimonious channel model such 
as ours. However, preliminary experimental results im, which 
use compressive measurements to successfully recover power- 
angle profiles for a controlled experiment (a small number 
of reflectors in an anechoic chamber), indicate that a simple 
ray-based model like ours may well suffice. Validating this 
assertion would require application of the more sophisticated 
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compressive estimation techniques discussed here, as opposed 
to the standard basis pursuit algorithms employed in Ifm . Note 
that existing models for simulation-based evaluations, such as 
the channel models standardized for IEEE 802.1 lad indoor 
60 GHz channels, typically assume more complex models 
which are variants of the Saleh-Valenzuela model, with a 
number of clusters, each consisting of multiple closely spaced 
rays. Such cluster-based models could be motivated by the 
roughness of reflecting surfaces such as walls, but they have 
not been experimentally demonstrated. While these issues fall 
beyond the scope of the present paper, which aims to make a 
fundamental contribution to signal processing for sparse spatial 
channels, a sustained effort in measurement-based validation 
of our model and approach is a critically important topic for 
future work. 

Alternative approaches to spatial channel estimation with 
RE beamforming include codebook-based techniques; in an 
indoor WPAN setting, it was shown in ifTSl that the minimum 
number of scanning beams for RE beamforming to get beam¬ 
forming gains close to (within a dB) that of exhaustive search 
is twice the number of antenna elements. This approach does 
not scale to the large arrays of interest to us. Hierarchical 
codebook search 03 is more efficient, but still requires far 
more overhead than our scheme, since it does not exploit 
channel sparsity. 

A hardware enhancement to pure RE beamforming as 
considered here (where a single RE chain serves all antenna 
elements) is to employ hybrid analog-digital beamforming, 
with a number of RE chains smaller than the number of 
antenna elements. We may term this an array of subarrays, 
with RE-level control for subarrays, and digital processing 
of subarray outputs. Such an approach is used in m, iisi 
for spatial channel estimation. Our work shows, however, 
that a single RE beamformed array suffices for this purpose. 
Of course, arrays of subarrays are certainly required for 
more advanced functionalities such as multiuser MIMO na, 
spatial multiplexing ([17], ([TSl, ([19], ||20|, IItTI . and spatial 
diversity im, I22|,|i23]. Integrating the compressive approach 
proposed in this paper within an array of subarrays architecture 
is an interesting area for future work. 

The present paper builds on our prior conference papers on 
compressive array adaptation HH, El, but goes well beyond 
them in several respects. In addition to a more detailed devel¬ 
opment of the analytical framework underlying our estimation 
algorithm, we now explicitly model the receive array at the 
mobile, which requires a generalization of the beaconing and 
feedback strategy. We also address system level design for 
compressive tracking in far greater detail, discussing link 
budget and overhead, and accounting for inter-cell beacon 
interference. Our initial work on compressive array adaptation 
II 24 I subsequently led to a general framework for compressive 
estimation Ea, m, which identify the isometries required 
to preserve fundamental bounds such as the Ziv-Zakai (ZZB) 
and Cramer-Rao (CRB), and use the relationship between 
these bounds to provide criteria for determining the minimum 
number of compressive measurements required to preserve 
geometry and to permit accurate parameter estimation based 
on a signal corrupted by an AWGN. We now adapt these 


general results here to develop guidelines for system-level 
parameter choices. 

HI. System Model 

Given the high demand for wireless data in dense urban 
environments, we focus our modeling and performance eval¬ 
uation on the urban canyon setting depicted in Eigure[T| with 
streets flanked by buildings on both sides. Picocellular base 
stations are deployed on lampposts in a zig-zag pattern on both 
sides of the street. We consider mm wave transmission on the 
downlink (for beaconing and downlink data) and LTE or WiEi 
at lower carrier frequencies on the uplink (for feedback and 
uplink data). In terms of channel estimation and tracking, this 
could be viewed as a worst-case assumption, since two-way 
transmission on the same mm wave band could potentially be 
exploited using channel reciprocity. Eor the east-west street 
shown, each base station has two faces, facing east and west, 
respectively. Each face can have multiple antenna arrays, but 
for simplicity, we consider a single antenna array for each 
face, used for both compressive beaconing and downlink 
data communication. Mobile stations are modeled as either 
pedestrians walking along sidewalks, or cars moving along 
the street. 

In our simulations, we model K = A dominant paths from 
base station to mobile in our simulations; the line of sight 
(LoS), and the single bounce reflections from the ground 
or the side walls. Some of these paths may be blocked by 
obstacles (diffraction around obstacles is limited for small 
carrier wavelengths). We ignore multiple bounces, since such 
paths get attenuated significantly, especially because each 
bounce sees a smaller angle of incidence than a typical single 
bounce. However, our compressive estimation algorithm does 
not use the preceding assumptions on number of dominant 
paths as prior information, and automatically discovers and 
tracks paths. 

Eor a regular d-spaced square Nm x Nuj antenna array 
and a point transmitter in the far-fleld, the channel seen by 
the array is a mixture of 2D sinusoids, each corresponding to 
a propagation path, and is given by 

l=K 

hn,,n = e C, 1 < m,n,<N,D, 

1=1 

where gi is the propagation gain of the (-th path, uji = 
are the spatial frequencies corresponding to the 
(-th path (w.l.o.g. we have assumed that the square array 
is placed in the x-z plane, with its sides aligned to the 
coordinate axes) and hm,n refers the channel seen by the 
m, n-th antenna element,. The spatial frequencies of the l- 
th path are given by u]x,i = 27r(d/A) s'mdi cost/); and u!z,i = 
27r((i/A) sind; sin^;, where d denotes the array spacing, A the 
carrier wavelength and {9i,(f>i) the inclination and azimuthal 
angles of the /-th path relative to x — z plane. We vectorize 
the 2D sinusoid ^gj{‘^xm+uj„n)^ 0 < m,n < Nid — l] and 
denote the resulting long vector by x(iVi£),u>), where 
ui = {ujx,ujz) is the frequency of the 2D sinusoid. Vectorizing 
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[hm,m 1 < m,n < Nm] in an identical manner gives us 

l^K 

h= '^gix{NiD,uJi). 

1=1 

Now, consider a base station transmitter with a regularly 
spaced 2D array of size Nt^m x and a mobile receiver 

with a regular 2D antenna array of size Nr^io x Nr^m- Let 
H denote the corresponding x channel matrix. 

Denoting by hi the ith row of this matrix, hf is the response of 
the receive antenna array to the ith transmit antenna. Denoting 
X {Nt^iD, by Xi {uj) and x (A^r-,in, <^) by x^ (uj), under the 
far-held assumption, it can be shown that 

l=K 

H = ^ gixt (ivf) (uj[), gi eC. (1) 

1=1 

Since we know the array geometries (in this case, a regu¬ 
larly spaced 2D array), an estimate of the x 

MIMO channel matrix H can be efficiently arrived at by 
estimating the spatial frequencies and the associated gains: 
{(gi,ujf,u;[), I = 1,..., K}. Such a parametric approach is 
far more efficient that direct estimation of individual entries 
of H, and enables us to drastically reduce the number of 
measurements required. 

IV. Compressive channel estimation 

We now describe our compressive approach for spatial 
channel estimation, which consists of a channel sounding 
strategy and an estimation algorithm which allows a base 
station to estimate the propagation gains {|pi|} and the spatial 
frequencies {w*} in parallel for a// mobiles in the picocell. 

A. Channel sounding 

The basestation sounds the channel using M compressive 
beacons. Each beacon is a known signal sent using a different 
set of transmit weights. The weights are chosen uniformly and 
independently at random from a small set of coarse phase 
shifts (for e.g. from the set {±l,±j}, where j = 1). 

We may therefore view each beacon as being transmitted 
from a different “virtual antenna” with a quasi-omnidirectional 
pattern. Each of the M transmit beacons are repeated L times 
by the basestation (see Eigure|2l). Eor each of these M transmit 
beacons, a mobile receiver employs L “virtual antennas” to 
measure the channel response, using receive weights chosen 
uniformly at random from {±l,±j}. Let y{m,n) denote 
the response at the (m, n)th receive element due to a given 
transmit beacon. Letting 6(m,n, fc) G {±l,±j} denote the 
weight for (m, n)th receive element for the fcth virtual receive 
antenna, the response seen by the fcth virtual receive antenna 
is given by 

r(k) = ^ b{m, n, k) x j/(m, n), 1 < k < L. (2) 

l<m,n<iVr,iD 

These measurements are used to construct the M x L Multiple 
Input Multiple Output (MIMO) “virtual channel” matrix V 
between the virtual transmit and receive antennas. Note that we 
do not require that the base station know the receive weights 


used by the mobile, or that the mobile know the transmit 
weights used by the base station. 

Denoting the vectorized version of weights of the i-th virtual 
transmit antenna by (a vector of shape x 1) and that 

of the j-th virtual receive antenna by bj x 1), the 

i,j-th element of V (the channel between the {i,j)-th virtual 
transmit-receive pair) is given by Vij = a^Hbj. Letting A = 
[ai ... a^]^ and B = [bi ... b/,]^, it is easy to see that 

V = AUB^. Using O, we have that 

l^K 

V = . (3) 

1=1 

The channel measurements are perturbed by i.i.d measurement 
noise, and are given by 

yi,j \/Zij ~ CAf(0,(J ), 

where Pe is the per-element transmit power. Letting Y and Z 
denote M x L matrices with their i,j-th entries given by yij 
and Zi j respectively, the “measured virtual channel” is given 
by 

l=K 

Y = s/Ky+Z = ^,Y,gi {Axt (a;*)) (Bx, {iol)f+Z. 


B. Feedback strategies 

Our goal is to track the mm-wave spatial channel as 
seen from the basestation, described by the parameters P = 
{(|pi| , uif) ,1 = 1,..., K}. To this end, every mobile in the 
picocell needs to feed back a portion of the measured virtual 
channel Y to the basestation. Lrom (0, we see that the 
information regarding the spatial channel as seen from the 
basestation, given by P, is available in the column space of 
Y. Building on this observation, we propose two feedback 
strategies: 

(i) The receiver feeds back the entire matrix Y. 

(ii) A more elaborate strategy involves feeding back Q < L 

strongest left singular vectors of Y scaled by their corre¬ 
sponding singular values, i.e., if Y = Yl\Zi with 

<Ti > 172 > ■ ■ ■ > cr^ > 0, the receivers feed back D = 
[ciUi • • • (TqUq]. This strategy identifies the Q-dimensional 
subspace of the column space of Y with maximum energy. 

V. Estimation Algorithm 

We now present an algorithm to estimate the parameters 
{{\gi \ ,ujf) ,l = 1,... ,K} characterizing the mm-wave chan¬ 
nel as seen from the basestation. The same algorithm applies 
for both forms of feedback discussed in Section IIV-BI the 
entire measured virtual MIMO matrix Y or the dominant 
weighted left singular vectors D. 

The fcth column of Y is given by 

l=K 

yk = +Zk, k = l,...,L (5) 

1=1 

where Zfc ^ CAf (O, denotes the fc-th column of Z and 

hi,k = sfPlgihl Xr We assume the weight sequence 
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{bfc} used to construct the receive virtual antennas at the 
receive antenna array is not available at the transmitter, and 
hence cannot jointly estimate and However, since 
{bfc, k = were picked in an i.i.d manner, we 

have that {hi^k, k = l,...,L} are i.i.d realizations of a 
random variable with ¥,\hi^k\^ = PeE\gih]^Xr {u>l)\ = 

PeN^iD \9i'^■ This allows us to estimate \gi'^ as follows: 

k—L 2 

P,\gi\^ = {l/{LNl^j,))Y,\hi,k ■ 

k=l 

From here on, in Section |V] we use the notation uji to refer to 
tij* and x{u}i) to refer to Axt(u)f). Thus, the measurements 
can be written as 

l=K 

yk = Y^ hi,k^ (uJi) + Zfc, k=l,...,L. (6) 

1=1 

We now provide an algorithm to estimate {uJi, {hi k}}- 


A. Single path 

We first present an algorithm for estimating a single path 
K = 1, which forms the basis for our sequential estimation 
algorithm for K > 1. Omitting the path index I in (|6]l, we 
have 

Vk = hkx{uj) + Zk, fc = l,...,L. 

Given that {z;} are independent realizations of CJ\f (O, cr 

the maximum likelihood (ML) estimates of u),{hk} are given 

by: 

k=L 

Cj,\hk \ = argminy^ ||yfc - /ifex(u>)||^ . (7) 

For any u), the optimal hk-s are given by least-squares 
estimates: 

(‘^) = (x(w),yfe)/||x(u;)f (8) 

where {x, y) denotes x^y. Plugging into (|7]i, the ML estimate 
of u) is given by: 

^ k=L 

w = argmax . |(x (u;), yfc)|^ (9) 

l|x(w)|| 

and the ML estimate of hk is given by h^. (u;). We employ a 
two-step algorithm to arrive the ML estimates: a “detection” 
phase followed by a “refinement” phase. 

Detection phase: Using M two-dimensional FFT computa¬ 
tions, we precompute x(a;) for frequencies of the form u) S 
$ = {(27rfc/T,27r(/T), 0<fc,(<T-l}, T = RNiD,t, 
where R is the oversampling factor. We pick the frequency 
from $ which maximizes (|9]l. The corresponding gains are 
given by hk = h^. (w). We remove the contribution of the 
newly detected sinusoid from the measured channel response 
and this residual measurement is given by 

Yk = yk-hkx{Q}). ( 10 ) 

(This residue is used for sequential detection for iT > 1, as 
discussed shortly.) 


Refinement phase: Our estimate from the detection phase is 
restricted to the discrete set <i> and consequently we do not 
expect uj to be equal to the ML estimate given by (|9]l (where 
the maximization is over [—27rd/A, 27rd/A]^ with d being the 
spacing between transmitter antennas). However, if we make 
the grid fine enough, the best estimate of w in $ is expected 
to be close enough to the optimal solution to allow refinement 
via local optimization. In order to do this, we first fix the 
gain estimates {hk} and refine only the estimate of the spatial 
frequency ui by seeking the minimizer of the ML cost function 


II ., 2 

^(u;) === ^ yfc - ^/cx(a;) 
k^l 

in the neighborhood of the current estimate Co using the 
Newton method. This involves evaluating the gradient vector 
Q (lj) and the Hessian matrix T-L{oj) of C{ijo) at the current 
estimate Co. The corresponding expressions are given by: 


= 

HijiLo) = 


dCjco) 

duji 

d^Cjco) 

dujidujj 


k=L 


= -2y^Jfil(rk,hk 


k=l 

k=L 


= - 2 ^] (rk,hk 


dx{Lo) 
doji 

d^x{Lo) 


hk 


k=l 

2 j dx(uj) dx(io) 
diOi ’ diOj 


" dcoidioj 


>, 1 <LJ <2 


where lo = [u;i,W 2 ]. {rfc}-s are the residual measurements 
given by (fTOl i. The Newton update for Co is 


Co-t^Co — H ^{Co)0{Co). (11) 

We follow this up by updating our estimates | (ifc | by plugging 
the new value of Co in ®, i.e., 

hk ^ hi (Co) = (x (w), yfc) j ||x (Co)\f (12) 

and modifying the residues ({r^}) accordingly using (fTOl i. The 
algorithm alternates between the updates in (fTTIi and (fTSli for 
a few iterations. 


B. Multiple paths 

We now build on the preceding single path algorithm 
for the general setting of iT > 1. Suppose that our cur¬ 
rent estimate of the sinusoids/paths is given by Vk = 
I {hi^k}'^ , ( = 1,..., iT|. The residual measurements 
corresponding to a set of estimated parameters V is given by: 

Vfe(T’) = yk- 'Eu,tAhi,kHv ■ 

Detect a new path: Assuming that the measurements y^ are 
given by the current residue Vk{Pq) (corresponding to the 
q detected paths), we use the single path algorithm in Sec- 
tion lV-Al to detect and refine a new sinusoid {hg+i^k})- 

Let Vq+i denote the new set of estimated parameters Vg U 
{(‘^g+i: {K+i,k})}- 

Refine existing paths: Once we add this new path, we refine 
the parameters of all q + 1 sinusoids in Vq+i one by one. 
Consider the parameters {Coi,{hi^k}) of the /-th sinusoid. 
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We use the refinement algorithm in Section IV-AI to refine 
{lji, {hi^k}) by assuming that the measurements are given 
by the residual measurements after excluding the sinusoid of 
interest, i.e, Vk{'Pg+i \ {(u);, {fi; fe})}). Sinusoids are refined 
in a round robin manner, and the process is repeated for a few 
rounds; 1 —>■ 2 (g + 1) —>■ 1 (g + 1). 

Stopping criterion: The algorithm continues to add newly 
detected paths as long as the payoff, determined by the total 
amount by which the residue decreases is above a threshold 
T. i.e., if 

k=L 

E {\Wk{'P,)f-\Wk{'P,+i)f) >r, (13) 

fc=i 

the (g + l)-th path is added and the algorithm proceeds by 
searching for a new path. On the contrary, if the reduction 
in total residue is smaller than r, the algorithm terminates 
and returns the prior estimate of the parameters Vg. For the 
stopping criterion, we use t = 30a^ log {20Nt^iD) in our 
simulations. This stopping criterion is empirically determined, 
such that the signal energy r is large enough to comfortably 
cross the ZZB threshold (so that we can expect the frequency 
estimate to be accurate). 

C. Tracking 

We sound the channel often enough so that between any 
two successive channel estimation cycles, the geometry of the 
mm-wave channel, given by the spatial frequencies {w;} of 
the paths, do not change “significantly,” even if the path gains 
{gi} do. This ensures that angle of departure estimates from 
the prior sounding round do not become stale over the course 
of the communication phase during which they are needed 
for beamforming purposes. For example, if we do not wish to 
tolerate a beamforming loss of 3dB or more, then our estimate 
from the previous round tb should be close enough to the 
current ui so that 

|(x(u)),x(w))|y||x(a;)f > 0.5 

over the entire communication phase. This condition is met if 
||u> — < 0.5 X {2tt/ Therefore, the estimates of 

spatial frequencies from the previous sounding round are good 
approximations of their current true value (within a DFT spac¬ 
ing of We exploit this by using {cbj, I = 1 ,... ,iT} 

from the prior round to initialize our algorithm (as opposed to 
using the empty set {}). We do this by constructing the matrix 
X = [x (cbi) ... X {<^k)] and setting hij to be the (z,g)-th 
entries of (X^X) ^ X^Y, where Y = [yi ... y^]. We 
refine all parameters in Vk = i = 1,. •. ,K} 

using the refinement algorithm in Section lV-BI before proceed¬ 
ing to seek for new paths using the algorithm in Section IV-BI 
Deleting weak paths: Paths estimated in prior rounds may not 
be viable at the current time instant (e.g, because of blockage). 
Therefore, we need means to remove such stale paths. We 
use the stopping criterion ( fT3l l as a means to delete weak 
paths. If deleting the path under question and optimizing other 
parameters increases the residue by an amount smaller than t, 
we delete the path permanently. Otherwise we keep the path. 


VI. Protocol Parameter Choices 

In this section, we give a principled approach to choosing 
parameters of the compressive channel estimation protocol, 
namely the number of unique transmit beacons M, the number 
of receive measurement weights L and the minimum effective 
Signal to Noise Ratio (SNR) needed for channel estimation, 
which we use to choose the sounding bandwidth Ws- We then 
turn to the question of how frequently the channel has to be 
sounded. In Section IVIIII we take two scenarios and apply this 
recipe to arrive at system level parameters for the protocol. 


A. Number of compressive transmit beacons 

Our goal is to estimate the spatial frequencies {ujf} from 
measurements of the form 


l=K 

Yk = ^ hi^kAyit (u;f) + z^, A: = 1,..., L. 

1^1 

The algorithm in Section |V] aims to estimate parameters 
by minimizing the ML cost function; 


7 = 1 


yj - Ei=f {<^i) 


E j=L 

i=i 


A X (w*) - /lijXt (w‘)) -b : 


where p denotes an estimate of parameter p. If the compressive 
measurement matrix A ensures that 


^ X ELf - hijxt (w|)) 


(14) 


M 




for relevant ({wi}, {d;i})-pairs, the cost structure of the es¬ 
timation problem is roughly preserved. Therefore, estimation 
using compressive measurements is similar to estimation with 
all N'^iu measurements (except for an SNR gain, given by 
M) iB^ . The main idea behind compressive sensing is that 
a “small” number of random projections can ensure that the 
geometry preservation condition in d is met with high 
probability. Such results build on the celebrated Johnson- 
Lindenstrauss Lemma (JL Lemma). 

JL lemma f2n\ : Consider a finite collection of points S C 
M" and an m x n matrix $ with its entries cho¬ 
sen in an i.i.d. manner from Uniform{±l} with m > 
(4 + 2/3) log |S'|/(e^/2 — e^/3), then with probability at least 


m(l 


e)< 


||(l)u - <i)v||^ 
l|u-vf 


< m(l -b e) Vu, V £ S. 


(15) 


It can be shown that this result extends to the setting 
where S C C" and the entries of 4> are chosen from 
Uniform{ ± 1, ±g }. 

Returning to the geometry preservation condition in (fT4l i. 
we see that when spatial frequencies are restricted to an over¬ 
sampled DFT grid G of size R = this condition 

reduces to a 2K-isometry property of the measurement matrix 
A relative to the basis X, where X is the x R matrix 

with its columns given by {xj (a;) ; a; € G}. For some fixed e, 
the matrix $ € enjoy the s-isometry property 
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M 

Fig. 3; Maximum and minimum values of 
||AXu||^y (M ||Xu||^) for different values of M, the 
number of transmitter beacons, across 5 x 10® random 
realizations of 8-sparse u. The basis X corresponds to the 
responses for a 32 x 32 array evaluated uniformly over a 
i? = 64 X 32^-sized grid 


for the basis B (of size p x n) if there exists a constant C > 0 
such that 

for all s-sparse u in C". It can be shown using the JL 
lemma (with arguments similar to those in ll28l l that if 
m = 0(se“^ log n), a randomly picked $ satishes this 
isometry property with high probability (w.h.p). Therefore, 
when the number of unique transmitter beacons scales as 
M = O (iTe“^logi?) = O log then the 2K- 

pairwise isometry criterion w.r.t the basis X is met by the 
randomly picked sounding matrix A (it can be shown that 
C = M for our choice of scale), thereby ensuring that 
the geometry of the spatial frequency estimation problem is 
preserved. 

We consider the example of the 32 x 32 transmit¬ 
ter array and plot the maximum and minimum values of 
(1/M) II 74Xu||^y/||Xu||^ from 5 X10® random realizations of 
a 2K = 8-sparse u in FigureOusing the 64-times oversampled 
DFT grid as the choice of basis X. We see that this ratio is 
within [—5,3] dB when M > 30. This illustrates that for 
estimating K = A paths using a 32 x 32 array, measuring 
the response corresponding to M = 30 random beacons 
approximates the effect of measuring all 32 x 32 = 1024 
antenna elements individually. 


B. Number of compressive receive measurements 

While we do not track {a;/, i = l,...,Ff}, the spatial 
frequencies at the receiver, we need to ensure that the set of 
measurements made at the receiver have sufficient informa¬ 
tion to estimate transmitter spatial frequencies. Suppose that 
||Bxr(u;/)|| Ki 0, it follows from hi j = (cu[) 

that all L measurements {y^, 1 < / < T} will have very small 
contributions from the i-th path, i.e., \hij\ ~ 0,1 < j < L. 



Fig. 4: Maximum SNR degradation 

Ta\n^ \\B:x.r{^)\\^!{L\\'X-r{<^)\\^) for the most favorable 
realization (from 10“^ runs) of an L x matrix B with 

Nr,ID = 4 

To see this observe that 

j=L j=L 

= Pe = Pe \9^f \\BXr (Oil 

i=i i=i 

( 16 ) 

When we restrict the receive spatial frequencies to an over¬ 
sampled DFT grid G of size R = it can be shown 

that for L = 0(\ogR) = O {log Nr, id), II-Bx^ (a;)||^ Ri 
L ||xr (u;)||^ = LN^ijj w.h.p (a direct application of the JL 
lemma (fT5T l with the set S being {x^ (tu) ; tu € G} U {0}, 
where 0 denotes the zero-vector). This ensures that 

j=L 

fv PeLN^^jj w.h.p, 

7 = 1 

thereby capturing power along the *-th path. We perform 
computations for the maximum power lost across spatial 
frequencies when using a 4 x 4 array and plot the results in 
Figure |4] This shows that around 5 carefully chosen projec¬ 
tions (we pick the best measurement matrix from 10"^ random 
instances) suffice to ensure that SNR degradation (relative to 
the nominal value of L) is no greater than 3dB for a 4 x 4 
receive array. 

C. SNR for successful estimation 

The preceding criteria delineate the regime in which the 
geometry of the estimation problem is preserved approxi¬ 
mately. We now turn to another factor which affects estimation 
performance, namely the SNR. Consider measurements of the 
form 

ym,n = 0<m,n< Nt,lD-l, ( 17 ) 

where Zm,n are i.i.d. CAf{0, a^) and spatial frequencies wi, u ;2 
and phase f are parameters to be estimated. The Cramer 
Rao Bound II 29 I (CRB) for estimating oj\ from measurements 
(fTTI i is given by C'(ct^) = 6/(SNR — l)), where 

SNR = ||xt(a;)||^ycr^ = (same expression holds 

for a; 2 )- Assuming an uniform prior over [0,27r)® for the 






Fig. 5: ZZB threshold SNR SNRth for estimating the frequency 
of a Nt,iD X Nt^iD sinusoid as a function of Nt^m 


parameters (wi, a; 2 , </>), the Ziv Zakai Bound (ZZB) with 
periodic distortion !^ for estimating wi evaluates to: 


ZiSMR) = l'Q(J 


SNR 1 - 


sin {Nt,iDh/2) 


7Vt_i£)Sin(/i/2) 


hdh. 


An indicator of the SNR needed for successful estimation 
is the convergence of the ZZB to the CRB Ehl . We use 
the SNR beyond which the ZZB is within O.ldB of the 
CRB as a measure of this convergence. We plot this ZZB 
threshold SNR for different values of Nt^io in Figure |5] e.g, 
SNRth = 16.04dB for an 8 x 8 array and SNRth = 16.13dB 
for a 32 X 32 array. 

The total energy Etot corresponding to the i-th path collected 
across the ML measurements {y^, 1 < j < T} is given by: 


3=L 

Etot —^ ^ ^ * 

j=i 

Using (flbl l in the above, we have that 

Etot = ||-4xt (u;‘)||^ WBxr (wDf Pe \gif 

= MLPNl^M\ 

where P = N^^j^Pe is the total transmit power supplied to the 
Nt.iD X Nt,iD antenna array. The above approximation holds 
when M and L satisfy the preceding geometry preservation 
criteria in Sections IVl-AI and IVl-B I respectively. The effective 
SNR of the i-th path is given by SNReff = Etot/cr^. It is 
important to note that the per-measurement noise variance 
is given by where is the noise variance 

per antenna element. Assuming no interference (which we 
account for in Section IVII-Db . = NoWs with Wg denoting 
the sounding bandwidth and Nq the thermal noise floor. The 
reason for the scale factor in the expression for cr^ is the 

following: Our measurement process consists of multiplying 
the received signal at each antenna (of which there are N'^uj) 
by phasors b{m,n,k) G {±l,±j} and adding the resultant 
signal (as per (|2|i). Since thermal noise seen by the 
isolated receive antennas are independent random variables, 
we have that tr^ = x (NqWs). Therefore, the effective 

SNR of the i-th sinusoid is given by: 

SNReff = MPP \g^\‘^/NoWs . 


This must exceed the ZZB threshold SNRth for successful 
estimation. Noting that ML/Wg is the time taken for channel 
sounding, the ZZB threshold SNRth gives us the means to 
evaluate the minimum overhead in time to estimate the channel 
for a given path gain \gi\^: 

Time taken = ML/Wg > SNRth Nq^P \gi'^ (18) 

The size of the picocell gives us a lower bound on \gif' and we 
later use this to guide us in choosing the sounding bandwidth 
Wg using (fTSI l. 

D. Sounding rate 

We round off the discussion on choice of protocol pa¬ 
rameters by giving a rule of thumb for the rate/frequency 
fs at which the spatial channel {u)i} needs to be re esti¬ 
mated. We use the estimated spatial frequency d> for beam¬ 
forming purposes in the time period between two channel 
sounding rounds (communication phase sandwiched between 
consecutive sounding phases; see Figure |7]i. Following the 
discussion in Section IV-Cl we have that if \\u>(t) — u! Iloo < 
Tt/Nt.iD throughout the communication phase, where u;(<) 
denotes the true spatial frequency and Cj the estimate from 
the prior sounding round, then the loss in SNR, given by 
||x((a;(i))||^y(xf (a;),Xf (u>(t))), is smaller than 3dB. If we 
assume that the closest user to the basestation array is at a 
distance R meters and that the maximum speed of a user 
in the picocell is given by Umax meters per second, then 
the maximum change (in terms of the £oo-norm) in spatial 
frequency Aw between consecutive sounding phases, spaced 
1 //b apart, is given by 2'itdvma^/fsR^- The worst-case 
geometry which achieves this bound is when the user is at 
a distance R along the bore-sight of the array and heading in 
a direction aligned with the one of the array axes. For this 
worst-case geometry (plotted in Figure |6]l, we have that: 

Aw < {2'Kd/X) sin A0 Ri tfUmax ^ 

JbXR 

Assuming that the estimate w from the previous sounding 
phase is accurate, if we ensure that 27rdumax//BPA < 
Tt/Nt^iD, we have that the beamforming losses in the inter¬ 
vening period are smaller than 3dB. This tells us that channel 
needs to be sounded often enough so that 

fB > 2dUmaxA't,lD/PA . (19) 

In the following discussions, we use the preceding in con¬ 
junction with (fTSl) to determine the overhead incurred in 
estimating the channel using the compressive architecture 
proposed herein. 


VII. System design 

We now discuss some key aspects of downlink system 
design related to our compressive architecture. We start by 
choosing basestation transmit power based on rules set by 
regulatory authorities, and then Ailing in the other details 
of the protocol according to the prescriptions laid out in 
Section |VT] Fixing the mobile array to be 4 x 4 {Nr, id = 4), 
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Fig. 6; Geometry corresponding to maximum change in 
The user moves along the 2 axis at a speed of Umax in the 
time interval 1 /Jb between two consecutive channel sounding 
rounds 


we consider two different choices for the base station array 
size; 8x8 (Nt^m = 8) and 32 x 32 = 32). All 

arrays are d = A/2-spaced. The total available bandwidth for 
communication and sounding is 2GHz. 

A. Transmit power 

We fix the effective isotropically radiated power (EIRP) 
to 40dBm, consistent with Federal Communications Com¬ 
mission (FCC) regulations for 60 GHz unlicensed transmis¬ 
sion. Accounting for transmit beamforming gain using an 
Nt.iD X Nt,iD array, the total transmit power 

P = 40 - 20 log Nt,iD dBm (20) 

which evaluates to 22dBm and lOdBm for Nt,iD = 8 and 
Nt.iD = 32 respectively. Assuming that this power is evenly 
split among the transmit elements, the power per trans¬ 

mit element is given hy Pe = P — 20 log W* id dBm = 40 — 
201ogA^t,iD dBm, which evaluates to 4dBm and —20dBm, 
respectively. Assuming that we design each element to operate 
at a fixed power, this is also the power per element used in 
the beaconing phase, even though the latter does not get the 
benefit of transmit beamforming. 

B. Communication range 

In order to ensure that the SNR for compressive estimation 
is adequate over a picocell, we first determine the picocell size 
using a nominal communication link budget, and then calculate 
the overhead required for successful estimation at that range. 
Standard link budget calculations, assuming oxygen absorption 
of 16 dB/km, an EIRP of 40dBm and a 4 x 4 receive array 
providing directivity gains of 12 dBi, can be used to show 
that we can attain a per-symbol SNR of 6dB at a link margin 
of lOdB for a symbol rate of 2GHz at a range of 200m. For 
omnidirectional free space propagation, the power gain in dB 
as a function of range r is given by 

GdB (r) = -p,r +20 logio 7 — (^B) 

47rr 

where /r = 0.016dB/m to account for oxygen absorption. 
Note that /i can be increased in order to account for rain. 
However, since our purpose is to ensure that channel estima¬ 
tion is successful whenever communication is successful, the 
contribution due to GdBij) cancels out, as we show shortly. 
Thus, while the particular value of p, determines picocell size. 


we shall see that it does not affect the overhead for channel 
estimation. 

The SNR per symbol is given by 
SNRc(dB) = EIRP (dBm) -f GdB(r) + 20 log^o 

’ ( 21 ) 

- lOlogio (NoWc) - Lraargin{cOmm) 

where Lmargin{comm) is the link margin (dB) for communi¬ 
cation. Note that lOlog^gWo = —174 -f NF dBm over a 
bandwidth of 1 Hz, where NF denotes the receiver noise 
figure in dB. Plugging in Wc = 2GHz, N^bd = 4, and 
NF = 6dB, we obtain a per symbol SNR of 7 dB at a range 
of r = 100 meters. 


C. Channel sounding protocol 

Our channel sounding protocol is specified by four param¬ 
eters: (i) bandwidth used by each basestation when sounding 
the channel, which we denote by Ws (ii) number of transmit 
beacons (or virtual transmit antennas) M (iii) number of 
receive measurements per transmit beacon (or virtual receive 
antennas) L and (iv) sounding rate /b which determines 
how often the channel is sounded. The parameters M, L 
and Wg together determine the effective sounding SNR. This 
must exceed the ZZB threshold SNR for successful channel 
estimation. This gives rise to the condition in [18] Imposing 
an estimation link margin Lmargin{est) (dB) and going to the 
dB domain, we have 


10 logio {ML/Ws ) > SNRth + ^margin {est) (22) 
-f 10 logio No- P- GdB{r) 

Adding (l2^ and (l2Tli and simplifying, we obtain 


10 logio (ML/Ws) > SNRth - SNRc -f L^argzniest) 

LmarginipOVri'ni) -\- 20 logiQ Ni iB 

-f 20 logio Nr, ID - 10 logio 


The key take-away is that ML/Wg must be large enough 
to compensate for the fact that we do not have the benefit 
of beamforming during the sounding phase. Notice that the 
range r (i.e., the dependence on picocell size) has cancelled 
out. Setting Lmargini^st) = 16dB (we use a higher link 
margin for channel sounding to account for power losses due 
to randomness of A and B), we obtain 


r.. . nil 

Time taken = - > 

Wg - 


16.34 ps 
0.2669 ms 


Nt,iD = 8 
Nt,iD = 32. 


(23) 


We choose the number of transmitter beacons for the 8x8 and 
32 X 32 transmitter arrays based on the geometry preservation 
criterion for the transmitter’s spatial channel estimation prob¬ 
lem discussed in Section lVI-AI We use M = 24 for Nt.io = 8 
and M = 30 for Nt,iD = 32 by numerically evaluating the 
worst-case distortion of pairwise distances relevant for the 
channel estimation problem (in Figure [3] we plot the worst- 
case distortion as a function of L for a random instance of 
A and Nt,iD = 32). Using the receive energy preservation 
criterion given in Section IVI-BI we choose the number of 
receive weights for the 4x4 receive array as L = 6. Using 
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Fig. 7: Channel sounding and communication phases of the 
proposed system 


these values for M and L in (l23T l. we obtain that the channel 
sounding bandwidth must satisfy 


< 


8.8124 MHz 
674.34 KHz 


Nt,iD = 8 
Nt,iD = 32. 


We choose = 8.8124 MHz for Nt^m = 8 and fF* = 
674.34 KHz for Nt^m = 32 so as to minimize the overhead 
in time, given by ML/Wg- Our specification of the channel 
sounding protocol will be complete when we give fs, the 
rate at which we sound the channel (see Figure |7]i which must 
satisfy (HSli. Assuming that the closest user is at a distance of 
R = 20m and that the maximum speed of a user in the picocell 
t^max is 45 miles per hour (20 m/s), we have that; /s > 8 Hz 
for Nt^iD = 8 and Jb > 32 Hz for Nt^io = 32. Choosing 
the minimum value for fs, we have that the overhead for 
our channel sounding protocol is MLfs/Ws = 0.0131% for 
Nt,iD = 8 and 0.8542% for Nt,iD = 32. 


allocated the same sounding BW) are well approximated by 
{kRfS, k £ 1i \ {0}}. Thus, the interference power seen by 
a single antenna is given by 

k—co k—oo 

^ = 2xEE PG{kRfS) = 8P E G(kRfS), 

k—1 paths k—1 

where we have assumed that there are 4 viable paths between 
the interfering basestation and user, each introducing the same 
amount of interference as the LoS path. This is a pessimistic 
assumption, since NLOS paths are attenuated by larger path 
lengths and reflection losses. Plugging in the expression for 
G(r), we have that 

k—oo 

I = /2t:'^R)S‘^) E 

k^l 

= {P\^/2t:‘^R)S‘^)U2 , 

where Li 2 (z) = jk^ is the dilogarithm function. 

The interference seen per antenna adds to thermal noise to 
give an effective per-element noise level of = NofFs + /. 
Assuming a worst-case geometry for the user of interest 
(distance of S from the basestation) and proceeding as in 
Section IVI-CI we see that effective Signal to Interference and 
Noise Ratio SINReg is given by 

SINReff= MLPG{S)/. 


D. Reuse analysis for channel sounding 

We investigate how a sequence of basestations employed 
in an urban canyon environment can share resources when 
estimating the spatial channel to users in their respective 
cells. The envisioned mm-wave system involves alternating 
between channel estimation and communication phases as 
shown in Figure [T] We assume that channel sounding rounds 
across basestations are aligned in time. We now characterize 
how the 2GHz spectrum is to be shared in space so as to 
limit the effect of interference from neighboring picocells on 
channel estimation performance. Such interference manage¬ 
ment is essential in the sounding phase; unlike the highly 
directive beams used in the communication phase, compres¬ 
sive sounding beacons are essentially omnidirectional. To see 
this, consider the average transmit power along any direction 
LO. This is given by Pe \\Axt{uj)\\'^^L « Pg ||xt(a;)||^ = 
PeNf^jj = P, the total transmit power. The approximation 
\\Ax.t{u})\\^!L K. ||xt(a;)||^ holds when the number of 
beacons L is large enough. Therefore, the average energy 
per-measurement received by an antenna at a distance r 
from a transmitter sending compressive beacons is given by 
PG{r), where G{r) = = A^/(l67r^r^) 

{u = (/r/10)lnl0) is the omnidirectional power gain at 
range r. We assume that basestations are deployed regularly 
as shown in Figure |9] and that the inter-basestation separation 
(along the street) is given by S. Suppose that the reuse factor 
is R (i.e, every Rth basestation uses the same slice of the 
frequency spectrum to estimate downlink spatial channels). 
We assume that for narrow urban canyons, the distance be¬ 
tween a user and all interfering basestations (those that are 


This can be rewritten as 

1/SINReff = 1/SNReff -b 1/SIReff , 

where SNReff = MLPG{S)/NqWs and the Signal 
to Interference Ratio SIReff = MLPG{S)/I = 
MLR)e-^^^. We need to ensure that 
SINReff exceeds the ZZB SNR threshold for successful es¬ 
timation. We choose the reuse factor Rf so that we are in the 
noise-limited regime by setting 

1/SIReff < 0.1 X (1/SNRth) « -10 - 16 dB. 

Assuming that protocol parameters are chosen so that SNReff 
exceeds SNRth, we can ignore interference in SINReff calcu¬ 
lations when 

SIReff > SNRth -b 10 Ri 26dB for Ntpo = 8, 32. 

In Figure |9] we plot achievable effective SIRs as a function 
of frequency reuse factor Rf for two example systems in 
Section IVII-CI i.e, 8x8 and 32 x 32 arrays with total 
number of measurements given hy ML = 24 x 6 and 
ML = 30 X 6 respectively. As the picocell size S grows, 
exponential attenuation due to oxygen absorption (the 
term in the expression for SIReff) helps in attenuating in¬ 
terference and improving SIR for same reuse factor Rf. To 
illustrate this we plot SIR as a function of Rf for three cell 
sizes S = 50,100, 200m in Figure|9] We observe that, in order 
to ensure SIReff > 26dB, a reuse factor of i?/ = 4 is needed 
for S = 50m, while Rf = 3 suffices for S = 200m. Plugging 
in the per-basestation sounding bandwidth Wg calculations in 
Section IVII-CI we see that the overall system-level channel 
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Fig. 8; Reuse of frequency resources for reuse factor i?/ = 3 



Fig. 9: Effective Signal to Interference Ratio SIReff for ML = 
24 X 6 (8 X 8 scenario) and ML = 30 x 6 (32 x 32 scenario; 
dashed) as a function of reuse factor Rf 

sounding bandwidth Wg 'x Rf for a S' = 50m picocell is as 
small as 8.8124MHz x 4 = 35.2MHz and 674.34KHz x 4 = 
2.7MHz for Nt^io = 8, 32, respectively, which is dwarfed by 
the total available bandwidth (2GHz). 

VIII. Simulation results 

We perform simulations for the two example systems con¬ 
sidered in Section IVIII (8x8 and 32 x 32 transmit arrays). 
We report results for the algorithm proposed in Section IV] and 
two feedback strategies: (i) ‘full’: users feedback the measured 
virtual channel matrix Y (M x L matrix; L = 6 for both 
systems) and (ii) ‘svd’: users feedback the 2 dominant left 
singular vectors of Y, scaled by their corresponding singular 
values (M x 2 matrix; one-third feedback overhead). 

We consider 6 mobile users moving in the urban canyon 
at speeds of 20, 3, 15, 1.5, 2.1 and 10 meters per second 
(covering both vehicular and pedestrian settings). The height 
at which each mobile device is held is in the 1.3 — 1.4m range. 
The basestation is mounted on a lamppost on the pavement 
(7 meters from a canyon wall), at a height of 6 meters. The 
basestation antenna array is tilted by about 7.5° in both the 
azimuth and elevation directions so that the boresight of the 
array points towards middle of the corresponding cell. This 
helps in more accurate spatial frequency estimation: since a 
change in direction near the boresight of the array results 
in larger changes in spatial frequencies than far away from 



Fig. 10: Six users in the urban canyon moving over the 
duration of the 7 second simulation interval. Their positions 
at time f = 0 is marked using a D-symbol 


the boresight, resolving paths is easier when the array points 
towards a direction in which we are likely to see more paths. 
We do not model blockage in these simulations, assuming 
that the LoS path and the three first order reflections are all 
available. Our goal is to estimate and track the K = 4 paths 
to all 6 users. 

Estimation error: Let T = {a>m : m = 1,..., K} denote 
the true spatial frequencies and P = {d>„ : n = l,...,K} 
denote the set of estimated spatial frequencies. When the 
base station uses one of the estimates in P, say d>, to 
form a beam, we do not realize the full 20 log iY* id dB 
beamforming gain. A measure of the sub-optimality is the 
estimation error ||u> — uj|| 2 , which we normalize by the DFT 
spacing of 27r/iVt id to define the following error metric: 

Aw(m) = min W^rn - '^n\\2 / (27r/iVt,iD) ■ (24) 

When no true spatial frequency exists near an estimate d>, i.e, 
when is a “phantom estimate”, we will quickly be able to 
discard it when we beamform in the direction of d> and And 
that the mobile does not receive power commensurate to what 
it expects with the 201ogiVt_iD beamforming gain. 

We plot the Complementary Cumulative Distribution Func¬ 
tion (CCDF) of estimation errors (EHi for the two systems 
{Nt,iD = 8, 32) in Figure [TTa] and the Probability Distribution 
Function (PDF) of the number of paths estimated K (correct 
value is K = 4) in Figure II Ibl From Figure [TT] we see that 
feedback of dominant singular vectors is an efficient feedback 
strategy which performs just as well as feeding back the entire 
matrix Y, while using only a third of uplink resources. 

Next, in order to evaluate the effect of errors in spatial fre¬ 
quency estimation on beamforming performance, we simulate 
a simple scenario in which the transmitter beamforms toward 
the strongest estimated path. Figure [12] shows the CDF of the 
achievable beamforming gain for an 8 x 8 array. While ideal 
beamforming requires adjustment of both gains and phases, 
suboptimal approaches for RF beamsteering with severely 
quantized phase-only control (four phases) have been studied 
in our earlier conference paper ll24l . We see from Figure [T2l 
that if ideal beamforming were performed with our estimates, 
then the SNR loss is less than 0.3 dB. If four-phase control 
is used based on our estimates, then the SNR loss is less than 
1 dB. The results for 32 x 32 arrays are entirely similar, and 
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(a) CCDF of frequency estimation errors dm for Nt, ID = S (left) 
and Nt,iD =32 (right) 


Fig. 11; Two feedback strategies considered: (i) M x 6 matrix 
overhead) 



Fig. 12: Beamforming gain achieved by an 8 x 8 antenna array 
for ideal and quantized beamforming techniques 


are therefore not plotted here.. 

Thus far, we have not said anything about channel frequency 
selectivity. Our proposed algorithm uses a small segment of 
the band to estimate the spatial channel, and the problem of 
channel dispersion is not addressed. However, we note that 
beamforming using a large array should reduce the effect 
of undesired paths, which simplifies the task of equalization. 
Figure [13] shows the channel impulse responses for the 32 x 32 
and 8x8 antenna arrays for a typical snapshot, when the 
transmitter beamforms towards the strongest estimated path. 
In our simulated setting, the LoS and ground reflection are 
close to each other in terms of both delays and angles of 
departure. We see that 8x8 array fails to resolve them, with 
both paths falling into the antenna’s main lobe, while the 
32 X 32 antenna array, which has smaller beamwidth (4° half 
power beamwidth), attenuates the undesired tap down to one- 
ninth of the desired path. Of course, it is possible to utilize 
the channel estimates far more intelligently, potentially with 
nulls directed both at strong undesired paths for the mobile of 
interest, and at the dominant paths for other nearby mobiles. 
The latter can be particularly useful for combating intra-cell 
interference when a base station face has multiple antenna 
arrays, each communicating with a different mobile. 



(b) PDF of # of estimated paths K (K = 4) for Nt,iD = 8 (left) 
and Nt,iD = 32 (right) 


Y (‘full’) and (ii) top two dominant singular vectors (one-third 



Fig. 13: Channel impulse response with quantized beamform¬ 
ing towards estimated strongest path for the 8x8 (left) and 
32 X 32 (right) scenarios 

IX. Conclusions 

We have shown that it is possible to super-resolve mm wave 
spatial channels with a relatively small number of compressive 
measurements, in a manner that is compatible with coarse 
phase-only control and RF beamforming. This allows scaling 
to a very large number of antenna elements without relying 
on channel reciprocity. While our discussion of system design 
issues such as link budget and inter-cell beacon interference 
is tailored to outdoor 60 GHz picocellular networks, the basic 
approach is broadly applicable (e.g., to other bands, and to 
indoor environments). An important topic for future work 
is comprehensive experimental validation of our compressive 
approach, which is based on a simple channel model including 
only the dominant rays. At the network level, there are a host 
of design issues (e.g., see discussion in 13). The compressive 
approach allows each base station to build up an inventory 
of viable paths to nearby mobiles, but there is a huge design 
space to be explored on how base stations coordinate using 
this information to alleviate the effects of blockage (mobiles 
in urban environments can be routinely blocked by pedestrians, 
automobiles, trees and other obstacles), and to manage inter- 
and intra-cell interference. Optimization of arrays of subarrays 
in base station “faces” for communicating with multiple users, 
as well as for handling channel dispersion, presents interesting 
design challenges. 
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